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Abstract. In the anomalous quantum Hall effect (QHE), a hallmark of graphene, nature of 
the edge states in magnetic fields poses an important question, since the edge and bulk should 
be intimately related in QHE. Here we have theoretically studied the edge states, focusing on 
the E = edge mode, which is unusual in that the mode is embedded right within the n = 
bulk Landau level, while usual QHE edge modes reside across adjacent Landau levels. Here we 
show that the n = Landau level, including the edge mode, has a wave function amplitude 
accumulated along zigzag edges whose width scales with the magnetic length, Ib- This contrasts 
with the usual QHE where the charge is depleted from the edge. The implications are: (i) 
The E — edge states in strong magnetic fields have a topological origin in the honeycomb 
lattice, so that they are outside the continuum ("massless Dirac") model, (ii) The edge-mode 
contribution decays only algebraically into the bulk, but this is "topologically" compensated 
by the bulk contribution, resulting in the accumulation over Is- (iii) The real space behavior 
obtained here should be observable in STM experiments. 



Ever since the anomalous quantum Hall effect (QHE) was experimentally observed, [TJ [2] 
fascination with graphene is mounting. The interests have been focused on the "massless Dirac" 
dispersions around Brillouin zone corners in graphene, where the Dirac cone is topologically 
protected due to the chiral symmetry [3]. The peculiar dispersion is responsible for the 
appearance of the n = Landau level (n: Landau index) precisely around energy E = 
in magnetic fields For the ordinary integer QHE an important question is how the bulk 
and edge QHE conductions are related for finite samples. Many authors have addressed this 
question[5j [6], and the bulk QHE conductivity, a topological quantity, is shown to coincide 
with the edge QHE conductivity, itself another topological quantity. This is an example of 
the phenomena that, when a bulk system has a topological order [TJ EJ EJ [10] that reflects the 
geometrical phase of the system [XT], this should be reflected and become visible in the edge 
states in a bounded system. [12| I13j This "bulk-edge correspondence" persists in graphene, as 
shown both from an analytic treatment of the topological numbers and numerical results for the 
honeycomb lattice [HI [To] . 

In this paper, we reveal the features in the real-space profile of the edges states in graphene 
in magnetic fields B in the one-body problem. A pecurier point on the graphene edge states is 
that the E = edge mode, despite being embedded right within the n = bulk Landau level in 
the energy spectrum, has a wave function whose charge is accumulated along zigzag edges. This 
is drastically different from the ordinary QHE, where edge modes reside, in energy, between 
adjacent Landau levels, and their charge is depleted toward an edge. The physics here indicates 
a topological origin in a honeycomb lattice, which is in fact totally outside continuum models. 
In B = 0, a zigzag edge in graphene has been known to have a flat dispersion at E = 0. [16] 
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Figure 1. Honeycomb lattice with armchair(a) or zigzag(b) edges with ei,e2 being respective 
unit translation vectors. Local charge density (oc area of each circle) in a magnetic field (ft = 1/41 
for the energy window —0.05 < E < 0.05 around n = Landau level is displayed, (c) Energy 
spectrum for a zigzag edge against hi (momentum along the edge) in a magnetic field of (ft = 1/41 
for zigzag edges, (d) Blowup around the n = Landau level for a zigzag edge, in a magnetic 
field of (ft = 1/5 here for clarity, with shaded regions represent the bulk energy spectra, while 
red curves the modes localized on the zigzag edge. 

which is protected by the bipartite symmetry of the honeycomb lattice. Here we focus on the 
edge states in strong magnetic fields, which has a flat dispersion at E = 0. 

We consider the tight-binding model on the honeycomb lattice with nearest-neighbor hopping 
t (which is taken to the the unit of energy hereafter). The magnetic field is introduced as a 
Peierls phase. The flux in units of the magnetic flux quantum is (ft = BSq/(2tt) = 1/q in 
each hexagon with an area Sq = (3\/3/2)a 2 . Since honeycomb is a non-Bravais, bipartite 
lattice with two sublattice sites • and o per unit cell, we can define two fermion operators 
c,(j) and c (j) with j = jiei + j2&2 defined in Figs. QJa) and(b) specifies the position of 
a unit cell. We assume that the spacing between the edges L\{j\ = 1,2, ■■■,Li) is taken 
to be large enough to avoid interference. The length along the direction (e2) parallel to the 
edge is also assumed to be long enough, for which we apply the periodic boundary condition. 
Performing a Fourier transform in that direction, we obtain a ^-dependent series of one- 
dimensional Hamiltonian, TL = J2k 2 ^iu(^2)- The resultant eigenvalue problem reduces to 
Hid(A;2)|V , (^2) E)) = E\ip(k2, E)), with corresponding eigenstates \ift(k2,E)). 

Having STM images in mind, we define the local charge density, 



Here x is the distance from the edge (as related to ji via ei which is not normal to the edge) , and 
Ei < E < E2 is the energy window to be included in the charge density (which is normalized 
to unity when the window covers the whole spectrum). 

Figure QJc) shows the energy spectrum for a zigzag edge in a magnetic field (ft = 1/41. For 
this magnetic field the n = Landau level around E = 0, with a narrow energy width, almost 
looks like a line spectrum on this plot. We calculate the local charge density defined in eq.([T]) for 
the armchair and zigzag edges. Figure [2{a) depicts the charge density I(x) normalized by the 
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bulk value Iq(= 4>) against the distance from the edge x measured by the magnetic length Ib for 
<p = 1/41, with the energy window \E\ < 0.05 set to cover the n = Landau level (along with the 
embedded E = edge mode). The charge density for an armchair edge decreases monotonically 
toward the edge, where the depletion occurs on the magnetic length scale {Ib = 3 3 / 4 a/y / 27r<^)), 
as in ordinary QHE systems. In sharp contrast, a zigzag edge has the charge density for the 
•-sublattice that is accumulated toward the edge while the charge density for the o-sublattice 
is depleted. When we perform this scaled plot the charge density with various magnetic fields, 
each of which coalesces on the common curves |17| . 

Outside the van Hove singurality, we recover the conventional Landau levels, so that we 
expect ordinary edge states. The result for the charge density for this energy region in FigJ^b) 
indeed shows that the charge density of outermost Landau levels (labelled B and C in FigUJ is 
depleted from the edge region over the magnetic length Ib in usual fashion. As expected, there 
is no difference in behavior between the A, B sublattices nor between armchair and zigzag edges. 

Going back to the n = Landau level for a zigzag edge, there is an important question: since 
the edge mode exists with an exactly E = flat dispersion, one might think that the charge 
accumulation along the edge entirely or primarily comes from this mode. We have checked this. 
We can use the transfer-matrix method [T4"] to examine the I{x) contributed by the E = flat 
band, for which nonzero amplitudes occur, rigirously, only on »-sublattice, for a sample with 
an infinite length along the edge direction. Figure Eta) presents I{x{ji)) for (f> = 1/21 and 
1/41, normalized by the bulk value of n = Landau level, Jo = 4>. Interestingly enough, we can 
immediately notice a plateau structure. In this plot where the horizontal axis j\ is normalized by 
g(oc 1/.B), where all the data points for different magnetic fields fall upon a common curve, which 
means that the chage density contributed by the edge mode has a series of plateau structures 
with a step arising every time ji increases by q for a magnetic field eft = 1/q. The height of the 
n-th plateau {I{x) with j\ = n) can be analytically given in terms of »-sublattice for the zero 
magnetic field, 

i ,1 +2(n-l) (i _ +2\ 

Asymptotically p n has an algebraic decay as p n ~ n~ 2 / {^r^/3). More interestingly, however, if we 
compare the total charge density with the contribution from the E = flat band in Fig. EKb) for 
the n = Landau level, the above plateau structure vanishes in the total density. This implies 
that, although the edge-mode contribution has a slow, algebraic decay, the bulk contribution 
compensates this, and we end up with the charge accumulation over the magnetic length scale. 
Since the Landau spectrum has a topological nature, we may call this curious phenomena a 
"topological compensation of charge densities", which is the final key result here. 

To summarize, we have shown that the charge density of n = Landau level in graphene in 
strong magnetic fields should be totally unlike ordinary QHE systems. The charge accumulation 
along zigzag edges only occurs for the E = edge mode in the n = Landau level, accompanied 
by a charge redistribution of the bulk states. The charge density around the edge can only be 
captured when bulk and edge contributions are considered simultaneously. The present result 
is expected to be measured by an STM imaging for graphene edges [15]. 
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Figure 2. (a) Scaled plot of the charge density I{x) for magntic flux <f> = 1/41 against x/Ib, the 
distance from the edge normalized by the magnetic length, for the n = Landau level (marked 
with A in Fig. [He)), (b) The same plot for the outermost Landau levels (marked with B and C 
in Fig. HKc)), f° r which Iq = 4>/2. 
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Figure 3. Scaled plot of I(x) at E = edge mode for </>(= 1/g) = 1/21 and 1/41 (a) and 
comparison with n = Landau Level contribution for <j> = 1/41. 
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